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Abstract. An elementary treatment of the Dirac Equation in the presence of a three- 
dimensional spherically symmetric S(r — r )-potential is presented. We show how to 
' handle the matching conditions in the configuration space, and discuss the occurrence 

of supercritical effects. 
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The problem of solving the Dirac equation in the presence of a 5 potential represents 
a curious and non simple situation, in contrast with the equivalent problem in non- 
relativistic quantum mechanics, i.e., the Schrodinger equation, which is discussed in any 
course on Quantum Mechanics. This is basically related to the fact that being the Dirac 
equation of first order, a singular potential, like the 5 one, induces discontinuities at the 
level of the wave function themselves instead of the usual discontinuities that appear in 
the first derivative in the Schrodinger scenario. 

This puzzling situation has been discussed previously in the literature by many 
authors. Rigorous construction of self adjoint extensions for the Dirac operator, allowing 
the handling of matching conditions at the support of the 5 potential were discussed in 
IJ. Using their results, in || a discussion of this problem was presented. However, we 
have realized that the proposed solution corresponds to a different self adjoint extension, 
associated to also a different singular potential. 

At present, a simple and elementary discussion, without the necessity of invoking 
sophisticated mathematical constructions, is still not available in the literature. In 
this note we want to overcome this situation, showing in elementary terms how to 
handle the problem in the configuration space. Other authors have shown how to handle 
the problem in momentum space 0. Using our results, we discuss the occurrence of 
supercritical effects, i.e., the possibility that the ground state starts to dive into the 
depths of the Dirac sea, implying positron emission ||]. 

The Dirac equation with an external potential can be written as 

Hip = Eip, (1) 

where 

H = ca ■ p + (3mc 2 + V(r). (2) 

Here a and (3 are the usual 4x4 Dirac matrices and ip is the four-component Dirac 
spinor. 

In what follows we will focus on the attractive spherically symmetric (vector) 
potential given by 

V(r) = —a8(r — r ), (3) 

with a > 0. 

It is convenient to recall some general properties of the solution of the Dirac 
equation in a central potential. For more details the reader may consult the book 
by Greiner, Miiller, and Rafelski || and reference therein. In this case, the complete 
set of commuting operators is given by H, J 2 , J3 and K, where J is the total angular 
momentum (i.e., J = L + S), and K is defined by 



K = /3(S- L + hl). 



(4) 



In this expression E 



a 
a 



where the ex's are the Pauli matrices. In terms of 



E, S = |E. On the other hand L is the orbital angular momentum. The eigenvalues of 
the operator K are given by 



where ip 



Kip 



(5) 



is the four-component spinor solution of the Dirac equation written 



in terms of the two-component upper and lower ip u and ipi. In the equation (|5]), 
j = 0, 1, |, ... are the usual eigenvalues of J 2 according to: 

jV = j(j + i)^V (6) 

We note that the four-component Dirac spinor is not an eigenfunction of L 2 . 
However the upper and lower components, taken separately, satisfy 

L 2 tfj u (x) = \j(j + l)h 2 + nh 2 + hi 2 ]^u{x) = l u {l u + l)h 2 4> u (x), (7a) 



and 
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L%(X) = + l)k 2 - Kk 2 + -?i 2 Wl(x) EE k(k + i)h%(x). 



(7b) 



Note that the orbital parities of the upper and lower components have opposite signs. 
It is convenient to parametrise the four-component spinor by separating the radial and 
angular dependence according to 







g(r)n jlum (6,(f)) 






_ if(r)£lji im (p,<p) _ 



where 



n 



jlm\ 



E 

m'm a 



(lsj\m'm s m) Y lm >(6, 4>)xi 



(8) 



(9) 



are the spherical spinors that carry out the angular part. Here we have coupled, through 
appropriate Clebsh-Gordan coefficients, the scalar spherical harmonics Y/ m /(6 l ,0) with 
the eigenfunctions of the spin given by the two-components spinors X\m 3 - 

For the present discussion we do not need the explicit form of the angular part. If we 
now consider the Dirac equation for a spinor parametrised in this way, it is not difficult 
to show that the radial components satisfy the following set of coupled differential 
equations: 



he ( — K — 

V dr r 



E -V(r) -mc 2 )G(r 



(10a) 
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and 



hc \^ + K j) = (E-V(r)+mc 2 )F(r), (10b) 



where G(r) = rg(r) and F(r) = rf(r). Because of the linear discontinuity of the spinor 
function, given by the delta potential, we need to fix the boundary conditions in the 
neighbourhood of the shell r = r . If we multiply equation (|1 qj) by F, equation (|10fl) 
by G, and then, sum both expressions to remove the singular delta function, we obtain: 

F ^ +G < G= ™ + K (£1^!> (11) 

n rnc 

In the previous expression, the primes denote radial derivatives. By integrating 
between r — e and r + e and taken then the limit when e — > we get: 

rr +e rr +e f 2mcFG ( F 2 - G 2 )\ 

lim / (FF + GG) dr = lim / I + k K „ ' dr (12) 

^ J ro-e e-^OJro-e \ % rTlC J 

Assuming that the discontinuities of these functions are well behaved we find that: 

lim(F 2 + G 2 ) |:^=0. (13) 

We may consider F and G as the real and imaginary parts of a function in (D . In 
this context, equation (|13|) expresses the fact that the absolute value of this function 
is constant when crossing the support of the 5-potential. This is in agreement with 
the condition established by J. Dittrich, P. Exner, and P. Seba (see equation (3.5a) in 

see also the remarks after equation (21) in 0). In fact, the absolute value of this 
function is continuous for all r. 

If we set = F(r ± e) and G + - = G(r ± e), (p~3f ) becomes 



F + 2 + G + 2 = FJ + GJ. (14) 

Now, as a second step, we multiply the differential equations by G and F, 
respectively, and substract them to get 

F'G - FG> = - (E T mC V + ( ^ + mC V + 2^1 - ^5 (r - r ) (f 2 + G 2 ) . (15) 
Tic he her he v ' 

Dividing by F 2 + G 2 , which is continuous for all values of r, we can integrate in the 

neighbourhood of the shell radius 

rro+e JP'Q _ JPQ> a rr +e 

lim / -— — — dr = — — lim / 5 (r — r ) dr. (16) 



e^O y ro _ E (F 2 + G 2 ) he e^O jr -e 

By using 

F'G-FG' _ 1 d /F\ 

[F 2 + G 2 ) ~ (F/G) 2 + lo> \G J ' (1 ' 



and since / 



l+h 2 (x) 



d[h(x)] = arctan(/i(a;)) we have: 
F(r)' 



lim arctan 



G{r) 



Iro— £ 



a 

he 



;is) 



In this way, our second boundary condition can be written as 



arctan 



arctan 



G. 



a 

he 



(19) 



Expressing the coupling constant in units of he, we define the dimensionless 
parameter a = tan(a/hc). Our second boundary condition can then be expressed as 



G 4 



{F-/G- 



a 



1 + a(F_/GL) 



(20) 



Except for an arbitrary phase, the last expression can be written as a matricial 
relation between the radial functions at both sides of the potential, 







cos(a/hc) 


— sm(a/hc) 




' F_ 






sm(a/hc) 


cos(a/hc) 




G - 



A 



(21) 



This matrix A is unitary (actually orthogonal), detA = 1, and contains the 
information for finding the eigenvalue equation for the bound states. Returning to our 
complex valued function with real and imaginary part given by F and G, respectively, 
it is interesting to remark that the 5 manifests itself through a change of phase of this 
function, given by t&n(a/hc). 

For the solutions of equations ( |10a|) and QlOfl) corresponding to the free case, we 
may separate the space into two regions: 

Region I, r < r , 



/(r)=^r(^-) J Jw+1/3 (At), (22) 
khc ( 7T y/2 



khc / 7T \ 

Fj(r) = Aj— -r — Ju.+Wfcr). (23) 



Region II, r > r , 



Gn(r) = A n r (^) 1/2 K lK+1/2 (kr), (24) 
F n (r) = -A n ™ C - r — K Ult+1/2 (kr). (25) 



E + mc 2 \2kr 
The relations we are looking for reduce to 
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F)_ khc Ii_ K+ i/ 2 (kr) 

Gi E + mc 2 h K+ x/2{kr) 1 ' 

Fn_ = _ khc K^ K+1/2 (kr) 

G n E + mc 2 K lK+l/2 (kr) 1 ' 

Taking into account that for the ground state, for j = I + s = 1/2, we have l K = 
and = 1 we can write 




In these equations k denotes a wave number, Tick = \/m 2 c 4 — E 2 . In order to find 
the eigenvalues of the hamiltonian, we evaluate (^6|) and (|27|) of r and use the boundary 
condition (|20|). In this way we are led to solving the following transcendental equation 



khc / 1 \ +a f khc Y (i 1 \ / 1 — tanh(fcr )A 
E + mc 2 \ kr J \E + mc 2 J \ kr J \ tanh(fcr ) / 

khc ( 1 — tanh(fcro) A 



E + mc 2 \ tanh(fcro) / 



a (32) 



Without solving explicitly this equation, we can analyse the behaviour of the ground 
state energy E, as a function of the parameter a = tan(a/frc), related to the coupling 
constant a. We start by introducing the following dimensionless variables 

E 



mc 



2 ■ 



(33 a) 



h/mc 

s = pvl — e 2 , (33 c) 

u = p(l + e), (33d) 
tanh(fcr ) 

90 = (33e) 
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In terms of these variables, our eigenvalue equation ( |32| ) can be written as 



s Q u Q (l + g s ) 
V - {so + l)s (l - #o) 

For a fixed value of the radius of the S -shell, ro (with r$ ^ 0), we are interested in 
determining the existence of the ground state energy in the interval (—mc 2 ,mc 2 ) (i.e., 
— 1 < e < 1). The existence of an e in this range will depend on the values of the 
coupling constant a through the dimensionless parameter a. 

Since Tick = ^/m 2 c A — E 2 = mc 2 \/l — e 2 , we get from (|33 e| ) that 

lim 5-0 = 1- (35) 

£— >±1 

On the other hand, using (|33dj we get 

lim - — — = \p. (36) 

Using the limits (|35"D and (^6]) in Eqn. (|3^) we see that as E approaches the free 
state (i.e., e — > 1 _ ), a — > + , which agrees with the fact that the bound states disappears 
for a vanishing potential. 

Proceeding as before, with e — > — 1 + (i.e., as the energy approaches the Dirac sea), 



we infer from equation (31) that a approaches the value — 3/2p, which gives 



where a crit , the minimum positive solution of (p7|) , is the value of the coupling constant 
for which the ground state energy sinks into the Dirac sea (e.g., for p — 1, a cri t = 2.19/ic). 

The numerical solution of equation (|34]) for the dimensionless ground state energy 
e as a function of p and a is plotted in Figure 1. Notice that for a fixed value of p, £ is a 
decreasing function of a. For all finite values of p there are supercritical effects. Clearly, 
the value of a cr u for which e sinks into the Dirac sea, is increasing with p. However, 
the limit r — > is not well defined in (^) and thus we cannot find solutions for bound 
sates in this limit. We would like to remark that a general theorem by Svendsen, || [7]], 
tells us that supercritical effects are absent in this limit. 
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Figure 1. In Fig.l, we show the behaviour of the ground state energy e = E/mc 2 , 
for different values of the size p — fo/(^) of the delta-shell as a function of the 
dimensionless coupling constant A = a/ (he): p = 0.5 (dot line), p — 1.0 (solid line), 
p = 2.0 (dashed line), p = 10.0 (dot-dashed line). 




Figure 1 



